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Jupiter and Saturn formed in a few Myr 1 from a gas-dominated protoplanetary disk and 



were susceptible to disk-driven migration on timescales of only ~100 Kyr. 2 Hydrodynamical 
simulations show that these giant planets can undergo a two-stage, inward-then-outward, 
migrationPS The terrestrial planets finished accreting much later, 6 and their characteristics, 
including Mars' small mass, are best reproduced starting from a planetesimal disk with an 
outer edge at ~1 AU. 7 8 Here we present simulations of the early Solar System that show how 
the inward migration of Jupiter to 1.5 AU, and its subsequent outward migration, leads to a 
planetesimal disk truncated at 1 AU, from which the terrestrial planets form over the next 
30-50 million years, with a correct Earth/Mars mass ratio. Scattering by Jupiter initially 
empties, but then repopulates the asteroid belt, with inner-belt bodies originating between 1- 
3 AU and outer belt bodies originating between and beyond the giant planets. This explains 
the significant compositional differences across the asteroid belt. The key aspect missing 
from previous models of terrestrial planet formation is an inward, and subsequent outward, 
migration of Jupiter. We conclude that the behaviour of our giant planets, characterized 
by substantial radial migration, is more similar to that inferred for extra-solar planets than 
previously thought. 

Hydrodynamic simulations show that isolated giant planets embedded in gaseous protoplane- 
tary disks carve annular gaps and migrate inward. 9 Saturn migrates faster than Jupiter; if Saturn is 
caught in the 2:3 mean motion resonance with Jupiter (see SI, Section 3, on conditions for this to 
happen), where their orbital period ratio is 3/2, generically the two planets start to migrate outwards 
until the disappearance of the diskP^ESl Jupiter could have migrated inward only before Saturn 
approached its final mass and was captured in resonance. The extents of the inward and outward 
migrations are unknown a priori due to uncertainties in the disk properties and relative timescales 
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for Jupiter and Saturn's growths. Thus we search for constraints on where Jupiter's migration may 
have reversed or "tacked". 

The terrestrial planets are best reproduced when the disk of planetesimals from which they 
form is truncated with an outer edge at 1 AU (refs 7,8). These conditions are created naturally if 
Jupiter tacked at ~ 1.5 AU. However, before concluding that Jupiter tacked at this distance a major 
question needs to be addressed: can the asteroid belt, between 2-3.2 AU, survive the passage of 
Jupiter? 

Volatile-poor asteroids (mostly S-types) predominate the inner asteroid belt, while volatile- 
rich asteroids (mostly C-types) predominate the outer belt. These two main classes of asteroids 
have partially-overlapping semimajor axis distributions,^^ though C-types outnumber S-types 
beyond ~2.8 AU. We ran a suite of dynamical simulations to investigate whether this giant planet 
migration scenario is consistent with the existence and structure of the asteroid belt. Due to the 
many unknowns about giant planet growth and early dynamical evolution we present a simple 
scenario reflecting one plausible history for the giant planets (illustrated in Fig 1). We provide 
an exploration of parameter space in the Supplementary Information that embraces a large range 
of possibilities and demonstrates the robustness of the results. In all simulations we maintain the 
fundamental assumption that Jupiter tacked at 1.5 AU. 

Figure 2 shows how the planet migration affects the small bodies. The disk interior to Jupiter is 
3.7 Earth-mass (M e ) equally distributed between embryos and planetesimals, while the planetesi- 
mal population exterior to Jupiter is partitioned among inter-planetary belts and a trans-Neptunian 
disk (8-13 AU). The planetesimals from the inner disk are considered "S-type" and those from the 
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outer regions "C-type". The computation of gas drag assumes 100 km planetesimals and uses a 
radial gas density profile taken directly from hydrodynamic simulations^ (see the Supplementary 
Information for a more detailed description). 

The inward migration of the giant planets shepherds much of the S-type material inward by 
resonant trapping, eccentricity excitation and gas drag. The mass of the disk inside 1 AU doubles, 
reaching ~2M e . This reshaped inner disk constitutes the initial conditions for terrestrial planet 
formation. However, a fraction of the inner disk (~14%) is scattered outward, ending beyond 
3 AU. During the subsequent outward migration of the giant planets, this scattered disk of S-type 
material is encountered again. Of this material, a small fraction (~0.5%) is scattered inward and 
left decoupled from Jupiter in the asteroid belt region as the planets migrate away. The giant planets 
then encounter the material in the Jupiter-Neptune formation region, some of which (~0.5%) is 
also scattered into the asteroid belt. Finally, the giant planets encounter the disk of material beyond 
Neptune (within 13 AU) of which only ~ 0.025% reaches a final orbit in the asteroid belt. When 
the giant planets have finished their migration, the asteroid belt population is in place, while the 
terrestrial planets require ~ 30 Myr to complete their accretion. 

The implanted asteroid belt is composed of two separate populations: the S-type bodies origi- 
nally from within 3.0 AU, and the C-types from between the giant planets and from 8.0-13.0 AU. 
The actual asteroid belt consists of more than just S- and C-type asteroids, but this diversity is 
expected to result from compositional gradients within each parent population (see Supplementary 
Information). There is a correlation between the initial and final locations of implanted asteroids 
(Fig. 3a). Thus, S-type objects dominate in the inner belt while C-type objects dominate in the 
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outer belt (Fig. 3b). Both types of asteroids share similar distributions in eccentricity and incli- 
nation (Fig. 3c,d). The asteroid belt is expected to have eccentricities and inclinations reshuffled 
during the so-called Late Heavy Bombardment (LHB)£^E]the final orbital distribution in our sim- 
ulations match the conditions required by LHB models. 

Given the overall efficiency of implantation of ~0.07%, our model yields ~1.3xlCT 3 M e of 
S-type asteroids at the time of the dissipation of the solar nebula. In the subsequent 4.5 Gyr, this 
population will be depleted by 50-90% during the LHB eventPEHand a further factor of ~2-3 by 
chaotic diffusion. 15 The current-day asteroid belt is estimated to have 6xl0~ 4 M e , of which 1/4 
is S-type and 3/4 is C-typeP^ Thus our result is consistent within a factor of a few with the S-type 
portion of the asteroid belt. 

The C-type share of the asteroid belt is determined by the total mass of planetesimals between 
the giant planets, and between 8-13 AU, which are not known a priori. Requiring that the mass 
of implanted C-type material is 3 times that of the S-type, and given the implantation efficiencies 
reported above, this implies ~0.8 M of material between the giant planets is left over from the 
giant planet accretion process, or ~16 M e of planetesimals from the 8.0-13 AU region, or some 
combination of the two. 

The simulations also found C-type material placed onto orbits crossing the still-forming ter- 
restrial planets. For every C-type planetesimal from beyond 8 AU that was implanted in the outer 
asteroid belt, 1 1-28 C-type planetesimals ended up on high-eccentricity orbits that enter the ter- 
restrial planet forming region (q < 1.0-1.5 AU; see Fig 3), and may represent a source of water 
for Earth. 16 For the Jupiter-Uranus region this ratio is 15-20 and for the Uranus -Neptune region 
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it is 8-15. Thus, depending on which region dominated the implantation of C-type asteroids, we 
expect that 3-1 1 x 10~ 2 M e of C-type material entered the terrestrial planet region. This exceeds 
by a factor of 6-22 the minimal mass required to bring the current amount of water to the Earth 
(~5 x 10~ 4 M e G-2), assuming that C-type planetesimals are 10% water by masP^ . 

Concerning the terrestrial planets, the migration of Jupiter creates a truncated inner disk match- 
ing initial conditions of previously successful simulations, 8 though there is a slight buildup of dy- 
namically excited embryos at 1 .0 AU. Thus, we ran simulations of the accretion of the surviving 
objects for 150 Myr. Earth and Venus grow within the 0.7-1 AU annulus, accreting most of the 
mass, while Mars is formed from embryos scattered out beyond the edge of the truncated disk. Our 
final planet mass vs. distance distribution quantitatively reproduces the large mass ratio between 
Earth and Mars and also matches quantitative metrics of orbital excitation (Fig. 4). 

Similar qualitatitive and quantitative results were found for a number of migration schemes, a 
range of migration and gas disk dissipation timescales, levels of gas density and planetesimal sizes 
(all described in the Supplementary Information). This scenario represents a paradigm shift in the 
understanding of the early evolution of the inner solar system. Here C-type asteroids form between 
and beyond the giant planets, closer to comets than to S-type asteroids. This can explain the vast 
physical differences between S-type and C-type asteroids, and the physical similarities between the 
latter and the comets (as shown by Stardust and micrometeorite samples see Supplementary 
Information for more on physical properties). 

If Jupiter and Saturn migrated substantially, then their birth region could be closer to the es- 
timated location of the snow line, ~3 AU, 21 rather than beyond 5 AU. Substantial migration also 
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points to similarities with observed extra-solar planetary systems for which migration is seemingly 
ubiquitous with giants commonly found at ~ 1.5 AUP^ However, the difference between our solar 
system and the currently known extra-solar systems is that Jupiter "tacked" at 1.5 AU to migrate 
outward due to the presence of Saturn. 
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Figure 1: The radial migration and mass growth imposed on the giant planets in the reference 
simulation. A fully-formed Jupiter starts at 3.5 AU, a location expected to be highly favorable for 
giant planet formation due to the presence of the so-called snow- lineP Saturn's 30 M© core is 
initially at ~4.5 AU and grows to 60 M ffi as Jupiter migrates inward, over 10 5 years. Inward type- 
I migration of planetary cores is inhibited in disks with a realistic cooling timescaleP^ thus, 
Saturn's core remains at 4.5 AU during this phase. Similarly, the cores of Uranus and Neptune 
begin at ~ 6 and 8 AU and grow from 5 M ffi , without migrating. Once Saturn reaches 60 M e its 
inward migration begins, 25 and is much faster than that of the fully grown Jupiter.^ Thus, upon 
catching Jupiter, Saturn is trapped in the 2:3 resonance. 3 Here this happens when Jupiter is at 1.5 
AU. The direction of migration is then reversed and the giant planets migrate outward together. In 
passing, they capture Uranus and Neptune in resonance which are then pushed outwards as well. 
Saturn, Uranus and Neptune reach their full mass at the end of the migration when Jupiter reaches 
5.4 AU. The migration rate decreases exponentially as the gas disk dissipates. The final orbital 
configuration of the giant planets is consistent with their current orbital configuration when their 
later dynamical evolution is considerecP^ (see SI section 3 for extended discussion). 
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Figure 2: The evolution of the small body populations during giant planets growth and migration 
as described in Fig. 1. Jupiter, Saturn, Uranus and Neptune are represented by large black filled 
circles with evident inward-then-outward migration, and growth of Saturn, Uranus and Neptune. 
S-type planetesimals are represented by red dots, initially located between 0.3-3.0 AU. Planetary 
embryos are represented by large open circles scaled by M 1 / 3 (but not in scale relative to the giant 
planets). The C-type planetesimals starting between the giant planets are light blue dots and the 
outer disk planetesimals are dark blue dots, initially between 8.0-13.0 AU. For all planetesimals, 
filled dots are used if they are inside the main asteroid belt and smaller open dots otherwise. The 
approximate boundaries of the main belt are drawn with dotted curves. The bottom panel combines 
the end state of the giant planet migration simulation (including only those planetesimals that finish 
asteroid belt) with the results of simulations of inner disk material (a < 2) evolved for 150 Myr 
(see Fig. 4), reproducing successful terrestrial planet simulations.^ 
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Distributions of 100 km planetesimals at the end of giant planet migration. Panel a: The semimajor 
axis distribution for the bodies of the inner disk that are implanted in the asteroid belt are plotted 
at three times: the beginning of the simulation (dotted histogram), at the end of inward planet 
migration (dashed) and at the end of outward migration (solid). There is a tendency for S-type 
planetesimals to be implanted near their original location. Thus, the outer edge of their final 
distribution is related to the original outer edge of the S-type disk, which in turn is related to the 
initial location of Jupiter. The final relative numbers of the S-type (red histogram, panel b), the 
inter-planet population (light blue) and the outer-disk (dark blue) planetesimals that are implanted 
in the asteroid belt are shown as a function of semi major axis. The orbital inclination (panel c) 
and eccentricity (panel d) are plotted as a function of semi major axis, with the same symbols used 
in Fig. 2. The dotted lines show the extent of the asteroid belt region for both inclination and 
eccentricity, and the dashed lines show the limits for perihelion less than 1.0 (left line) and 1.5 
(right line). Most of the outer-disk material on planet-crossing orbits have high eccentricity, while 
many of the objects from between the giant planets were scattered earlier and therefore damped to 
lower eccentricity planet-crossing orbits. 
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Figure 4: Summary of the 8 terrestrial planet simulations showing the mass vs. semi major axis of 
the synthetic planets (open symbols), compared to the real planets (filled squares). The triangles 
refer to simulations starting with 40 embryos of ~ 0.051 M®, and squares to simulations from 
80 embryos of ~ 0.026 M e . The horizontal error bars show the perihelion- aphelion excursion 
of each planet along their orbits. The initial planetesimal disk had an inner edge at 0.7 AU to 
replicate previous workj^and an outer edge at ~ 1.0 AU due to the truncation caused by the inward 
and outward migration of the giant planets as described in the text. Half of the original mass of the 
disk interior of Jupiter (1.85M©) was in ~727 planetesimals. At the end of giant planet migration 
the evolution of all objects inward of 2 AU was continued for 150 Myr, still accounting for the 
influence from Jupiter and Saturn. Collisions of embryos with each other and with planetesimals 
were assumed fully accretional. For this set of 8 simulations the average Normalized Angular 
Momentum DeficiP was 0.001 1 ±0.0006, as compared to 0.0018 for the current solar system. 
Similarly, the Radial Mass Concentration 30 was 83.8±12.8 as compared to 89.9 for the current 
solar system. 
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1 Outline 

Here we provide more details on the simulations and on the dependence of the results on the 
parameters of the problem. We discuss, in sequence, the following points: 

• Section 2: Simulation methods, for reproducibility of the results. This includes: (i) 
adopted gas profile, (ii) formulae for gas-drag and tidal damping calculations and (iii) 
implementation of giant planet migration. 

• Section 3: Defining a plausible evolution scenario for the giant planets. 

• Section 4: Relevance of different aspects of giant planet evolution for sculpting the inner 
solar system. 

• Section 5: Dependence of the results on the initial location of Jupiter and different migra- 
tion timescales. 

• Section 6: Dependence of the results on the size of the planetesimals. 

• Section 7: Implications for the Asteroid Belt. 

• Section 8: Placing this result in context. 

These issues underline the robustness of our results, but also its parametric dependencies. 

Our present paper should be considered as a conceptual demonstration that the inward- 
then-outward migration of Jupiter with a tack at 1.5 AU can explain properties of the inner 
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solar system that were unsolved problems until now. We did not attempt any quantitative fine- 
tuning of our results, but limited ourselves to showing qualitative agreements and appropriate 
efficiencies. The parametric dependencies of the results illustrated in Sections 5 & 6, however, 
should be useful guidelines for future work, aimed at a better quantitative reconstruction of the 
asteroid belt structure. 

Section 7 provides a discussion on the possible implications this work has related to the 
origin of the asteroid belt, and notes recent related research. We end this paper in Section 8, 
with some considerations to clarify how the present results should be viewed in the context of 
the evolution of the solar system after the removal of the gas disk, during the post-gas-depletion 
phases. 

2 Methods 

Gas Profile A specific gas profile is used, taken from a hydrodynamical simulation 31 of the 
evolution of Jupiter and Saturn in a gas disk with an aspect ratio of 5%. In that simulation, 
the gas surface density was 100 g cm 2 at 15 AU (in between two common estimated values 
32,33) anc ^ i nwarc [ s of this radius, it was carved by the presence of the two planets, with Jupiter 
at 5.2 AU (solid curve in Fig. SI). However, we are modeling the migration of Jupiter from 
its starting location down to 1.5 AU and then out to ~5.4 AU. Thus, the gas profile needs to 
change as the planets migrate. We use the fact that the gas profile must be invariant with the 
location of Jupiter if the latter is assumed to be the unit of distance; moreover, assuming that 
the unperturbed surface density of the disk is proportional to 1/r, we scale the surface density 
profile by the factor 5.2/r Jup , where r Jup is the orbital radius of Jupiter in AU. This gives, for 
instance, the dashed and dotted curves shown in Fig. SI when Jupiter is at 1.5 AU and 3.5 AU 
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respectively. 

There is no consensus on what the surface density of the disk was when Jupiter and Saturn 
formed. It might have been a few times larger or smaller than the Hayashi-Weidenschilling 
values. For the dependence of the results on the gas density, we invite the reader to check in 
Sect. 6 where we discuss the dependence of the results on the size of the planetesimals. In fact, 
under the effect of gas-drag, the evolution of a planetesimal of a given size in a disk /-times 
more massive is the same as that of a planetesimal /-times smaller in the nominal disk. Thus 
by testing multiple planetesimal sizes we also explored a similar range in gas density. 

During the outward migration of the planets, the gas profile illustrated in Fig. SI is mul- 
tiplied by a factor exp(— t/r), where t is the time from the beginning of the outward migration, 
and r is the timescale of the migration. This mimics the disappearance of the disk as the planets 
are moving, which ultimately brings the planets asymptotically to their final resting orbits. 

Gas Drag and Tidal Damping For gas drag purposes, we focus on 3 different fiducial plan- 
etesimal sizes: 10, 100 and 1000 km in diameter, each in a different simulation. We do not 
simulate smaller planetesimals, because current asteroids of 1 km in size or smaller are proba- 
bly not primordial, but rather are fragments of smaller objects during the subsequent collisional 
evolution 34 . 

The bulk density of the planetesimals is assumed to be 1.5 g cm -1 . Given their sizes, 
this gives their inertial mass for gas-drag calculation. Notice that the dynamical mass of the 
planetesimals in the S-type disk (used to compute their gravitational effects on the planetary 
embryos, which eventually will form the terrestrial planets) is in all cases 0.00038743M ffi . In 
this way, given the total number of 4830 planetesimal "tracers" in our simulations, the total 
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planetesimal mass in the S-type disk is 1.87M®. Conversely, the dynamical mass of the plan- 
etesimals from the C-type disk is set more than two orders of magnitude smaller, 1.3833e-6 

M e . 

We adopt commonly used formulae for aerodynamic gas drag 35 , and their implementation 
into the numerical integration code has been previously detailed 36 . To compute the volume 
density of the gas and its vertical profile, we assume that the scale height of the disk is 3.3% at 
1 AU and flares out as r 25 . 

Planetary embryos in the S-type disk also undergo tidal damping 36 . These formulae also 
depend on the local aspect ratio of the disk, which we assume as indicated above. The effect of 
type I migration on the evolution of the semi major axes of the embryos was not included, given 
the recent uncertainties on the real behavior of type I migration in the inner part of realistic 
disks 37 (see Section 3). 

Migration of the planets Each planet is forced to migrate by imposing a small change Sv to 
the orbital velocity along the velocity vector at each time-step. In this way, the norm of the 
velocity vector evolves with time. 

We choose two functional forms for the velocity change for inward and outward migra- 
tion. In the nominal simulation presented in the main paper they are: 

• for the inward migration 

v(t) = v + Av in [t/t CQd ] 

where v is the original velocity, t cnd is the time at which inward migration ends (10 5 yr 
in the simulation presented in the main text) and t is counted from the beginning of the 
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simulation; 
• for the outward migration 

v(t) = v- m + Av out [l - exp(-t/r)] 

where v- m is the velocity at the end of the inward migration (i.e. v + Ai> in ), r is the 
adopted migration timescale (10 5 yr in the simulation presented in the main text, but 
different values are tested below) and t is measured from the beginning of the outward 
migration. 

The values of Av - m and Av out are set from the locations of the planets at the beginning 
and at the end of each piece of migration (resp. inwards and outwards). The small velocity 
kick 5v , that is given to each planet at each time-step during the integration, is simply computed 
from the formulae above by differentiation relative to time. 

Because the inward migration of Jupiter might slow down significantly before the out- 
ward migration starts, in Sections 5 & 6 we adopt an exponential functional form for the ve- 
locity change during inward migration, analog to that used for outward migration. Whenever 
migration is exponential, the total integration time is 5t. 

3 Defining a plausible evolution scenario for the giant planets 

The accretion of giant planets is still one of the most mysterious phases in the history of the solar 
system and of planetary systems in general. There are, nevertheless, several known aspects of 
the process of planet formation and dynamical evolution that can help us to envision a plausible 
coherent scenario, like the one proposed in Fig. 1 of the main text. We review these aspects 
here. 
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Until a few years ago, it was thought that planetary cores had to migrate very fast towards 
the central star, as a result of their gravitational interaction with the gaseous component of the 
proto-planetary disk in which they form: the so-called Type I migration 38 ^ 2 . Type-I migration 
was a serious problem for the classical core-accretion model of giant planet formation. 

However, Type-I migration is valid only in isothermal disks. Instead, it has recently been 
shown that migration can be outwards in disks that transport and dissipate heat inefficiently 43 ^* 5 . 
This is an authentic paradigm shift in our understanding of how giant planets can form. Outward 
migration of planetary cores cannot proceed indefinitely, as the outer part of the disk can radiate 
heat efficiently, approaching the properties of an isothermal disk. Recent work has shown that 
there are two equilibrium radii in the disk where core-migration stops 46 . The innermost one 
of these equilibrium positions is associated with the snow-line (the location where water- vapor 
condensates in ice), supposedly around 3-5 AU. The location where core-migration stops has 
been shown to be a privileged site for the growth of cores 47 ' 48 . Thus, cores are expected to 
grow in the vicinity of an equilibrium radius, on non-migrating orbits, locked with each other 
in mutual resonances. This is what we assume in the scenario illustrated in Fig. 1 of the 
main paper, where the planets form in the 3-8 AU region and don't migrate until they reach a 
sufficiently large mass (about half a Saturn mass). 

Once planets achieve a mass of the order of half a Saturn mass, normal Type-I migration 
resumes 45 . These planets will then leave the equilibrium radius where they grew (see also ref 
46). If the disk is sufficiently massive, inward migration for Saturn-mass planets can turn into 
Type-Ill migration, which is a much faster, accelerating and self-sustained migration mode 49 . 
More massive planets like Jupiter, instead, open wide and deep gaps in the proto-planetary disk, 
and their migration enters the Type-II mode 50 - 51 . Type-II migration proceeds at the velocity at 
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which the gas is driven into the star by its own viscosity 52 and therefore it is, in general, much 
slower than Type-I or Type-Hi migration. 

Hydro-dynamical simulations, where Jupiter and Saturn are simultaneously taken into 
account with fixed masses, show that Saturn has no trouble in approaching Jupiter in a few 
hundred orbital periods and leads to the two planets being caught in resonance 53 . In disks with 
mass comparable to the so-called Minimal Mass Solar Nebula 32 ' 33 , the most likely end-state is 
the capture in the mutual 2:3 mean motion resonance, where the orbital period of Saturn is 1.5 
that of Jupiter. This occurs even if Saturn is initially beyond the 1:2 resonance or locked into 
the 1 :2 resonance 54 . Stable capture into the 1:2 resonance is possible only for disks with surface 
density decaying less steeply than 1/r or in low-mass disks 55 . 

Once locked in the 2:3 resonance, the inward migration of Jupiter and Saturn stops; the 
typical evolution is now outward migration 53 . The speed of outward migration depends on the 
depth of the partial gap opened by Saturn, which in turn depends on the properties of the disk 
(scale height and viscosity) 31 . For thick disks (about 6% in scale height for a typical viscosity), 
the migration is very slow, and the planets remain at effectively constant distance from the 
central star. But for disks with decreasing thickness, outward migration becomes increasingly 
fast. In principle, outward migration can bring the planets ten times further than their initial 
location in a few thousands orbital periods 56 . The peculiar dynamics of Jupiter and Saturn in 
the 2:3 resonance was immediately identified as the key to explain why, in our solar system, 
Type-II migration did not bring Jupiter closer to the Sun in stark contrast to most of the extra- 
solar giant planets detected so far 53 . Until now, outward migration of Jupiter and Saturn from 
inside ~4 AU was considered incompatible with the existence of the asteroid belt - therefore 
only proto-Solar disk models which prohibited inward migration were considered viable 31 ' 57 . 



7 



As we show in the main paper and in the next sections, however, this is not true, releasing the 
constraints against Jupiter's outward migration. This opens a new degree of freedom to model 
the evolution of the inner solar system. 

It is worth stressing that the phenomenon of migration reversal of two planets in reso- 
nance occurs only if the outer planet is less massive than the inner planet 31 , with the ideal mass 
ratio being between 1/4 and 1/2. The resonance that can drive outward migration is the 2:3 for 
Jupiter-Saturn mass planets and the 1:2 for more massive planets (for instance 3 times Jupiter 
and Saturn-mass, respectively for the inner and outer bodies). Instead, if the planets have com- 
parable mass, or the outer planet is more massive, the resonant planets migrate inwards 58 60 . 

Thus, all the ingredients to define the inward-then-outward migration scenario illustrated 
in Fig. 1 of the main paper seem to be established: (i) Jupiter migrates inward while it is alone 
in the disk (i.e. when Saturn is still a core); (ii) when Saturn grows towards its final mass it 
approaches Jupiter until it is trapped in resonance; (iii) when the two planets are in resonance, 
they migrate outwards. 

However, things are not so simple. All the hydrodynamical simulations that show that 
Saturn captures Jupiter in resonance, assume fixed masses for the planets. A natural question 
arises: is it still reasonable to expect resonance capture if the migration histories of the planets 
are coupled with their accretion histories? 

At first sight the answer is negative. If the second planet forms later than the first one, 
it's migration history should just replicate that of the first planet but later in time. More simply, 
the second planet should always lag behind the first one, as it is just repeating the evolution of 
the first planet, just at a later time. Thus, it appears that Saturn could catch Jupiter in resonance 
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only if the accretion histories of the two planets were different. In particular, if Jupiter grew 
very rapidly from half a Saturn mass to close to its current mass, it would have passed directly 
from no-migration to Type-II migration, which is relatively slow. To some extent this must 
have happened because Jupiter could not have spent considerable time at a Saturn mass or else 
it would have migrated very close to the Sun (due to the formidable migration speeds at this 
mass), which is not what we observe today. Instead, if Saturn grew more slowly than Jupiter 
and spent longer time near a Saturn-mass, it would have undergone Type-I or Type-Ill migration 
and hence would have approached Jupiter until being trapped in resonance. So, the question 
becomes whether it is reasonable to expect different accretion histories for the two planets. 

One could naively think that Saturn grew slower than Jupiter because there was less gas 
in the system at the time when Saturn formed. This is not a valid answer. In fact in first 
approximation both accretion and migration are linear in the surface density of the disk S, so 
the two effects cancel out: the amount of radial migration in the time interval corresponding to 
two given masses of the planet would be the same whatever the mass in the disk. In practice, 
the time t needs to be renormalized as t' = We need to find a mechanism for which Saturn 
grows slower than Jupiter in normalized time t'. 

Unfortunately, the accretion of gas by a forming giant planet is not a well understood 
process. Thus, it is unclear how accretion speed relates to several parameters of the disk, 
which might have changed in the interim between Jupiter's accretion and Saturn's accretion. 
There is, however, one well recognized parameter, that greatly impacts accretion speed: the 
disk opacity. The disk opacity is likely to have increased as time passes, because H and He 
are photo-evaporated and removed from the disk, while heavier molecules are not removed, or 
are removed at a slower rate given their larger masses 61 - 62 . This process has been invoked to 
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explain why Jupiter is enriched in elements heavier than Helium by a factor of 3-4 relative to 
solar composition 62 63 ; this assumed that Hydrogen and Helium had already been depleted by 
a factor of 3-4 in the disk by the time Jupiter captured its atmosphere. Saturn's atmosphere 
is enriched by a factor of 11 in Carbon 64 . Therefore, if this interpretation of the enrichment 
of the giant planet atmospheres is correct, one can conclude that the disk was richer in heavy 
elements when Saturn formed, and hence it had a larger opacity. The higher the opacity, the 
slower the capture of gas in the extended atmosphere of the planet 65 66 . Thus, this may be a 
valid explanation, although probably not unique, of why Saturn's gas accretion was slower than 
Jupiter's. 

Admittedly, given the unknowns on giant planet accretion, we cannot claim that Saturn 
must have captured Jupiter in resonance (given the great diversity of planetary systems discov- 
ered so far, the verb "must" probably does not apply to planetary evolutions). But the results 
and the arguments reported above definitely argue that it is plausible that Saturn did so. 

We now turn to extra-solar planetary systems for evidence about the possible evolution of 
a pair of planets. There are a few systems (for instance HD 134987) with a more massive inner 
planet and a less massive outer planet with a large separation and with the inner planet less than 1 
AU from the central star. In these systems the outer planet, despite its smaller mass, obviously 
never caught the inner planet in resonance. But the final orbital structure of these systems 
(large orbital separation; warm Jupiter) is very different from our own. So, it is evidence that 
something else had to happen in our system (or in the system around OGLE-06-109L, which is 
a "twin" of the Jupiter-Saturn system). There are many pairs of planets in resonance (or close 
to), near their star; none of these cases exhibits a Jupiter/Saturn mass ratio. The absence of 
this configuration, which is statistically significant if one assumes that the mass ratio should 
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be random, strongly supports the theoretical result that resonant planets in Jupiter/Saturn mass 
ratio move outwards, and therefore cannot be found within the range of stellar distances that 
can be probed by radial velocity observations (OGLE-06-109L's system was in fact discovered 
by micro-lensing). Finally, the HD 45364 system is probably in 2:3 resonance 67 . Fast inward 
migration of the outer planet is needed to overcome the 2:1 resonance 68 . If both planets had 
undergone this fast migration, there would be no planets left in this system. This indicates 
that planetary systems can survive with only one planet undergoing fast migration. However, 
the mass ratio in the HD 45364 is 1/3, so the analogy with the Jupiter- Saturn system cannot be 
pushed too far (and indeed these planets are close to their star, so they did not migrate outwards). 

In summary, the extra-solar planets cannot give us any definitive answer, but it is likely 
that there is a great source of diversity for the fate of Jupiter/Saturn-like planets: when the 
second planet does not capture the first one in resonance, the system ends up like HD 134987; 
if the opposite occurs, the system ends up like our solar system. 

To conclude this section, we note that the outward migration of Jupiter and Saturn is 
probably the most natural explanation for the compact, resonant configuration of the four giant 
planets (including Uranus and Neptune, see Fig. 1 of the main paper) that is assumed by the 
Nice model at the time of disappearance of gas 69,70 . The Nice model is extremely successful 
in reconstructing the current structure of the solar system, starting from dynamically cold and 
compact orbits (see last section). The multi -resonant configuration of the giant planets was 
a natural outcome of the evolution of the planets in the gas-disk because that work assumed 
that Uranus and Neptune were affected by Type-I migration 31 . As explained above, our current 
understanding has changed, and Uranus and Neptune should not have experienced migration. 
Thus, if Jupiter and Saturn had migrated inwards, without migration reversal, while Uranus and 
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Neptune did not migrate, at the disappearance of the gas the system should have had a large 
orbital separation between Saturn and Uranus/Neptune. Now, other works have excluded this 
configuration 71 , by showing that it cannot reproduce the current orbital architecture of the giant 
planets. 

In summary, migration reversal is required to explain, together with the Nice model, the 
structure of the outer solar system; moreover, as we show in this paper, it is required to explain 
the structure of the inner solar system. This unprecedented level of success argues strongly that 
resonant capture and migration reversal are real features of the true evolution histories of Jupiter 
and Saturn. 

4 Relevance of different aspects of giant planet evolution for sculpting the inner solar 

system 

As discussed above, there are large uncertainties for how the planets accreted and dynamically 
evolved in the gas disk. The scenario presented in the main paper and justified in the previous 
section is definitely not unique, even among those leading to the required tack of Jupiter at 
1.5 AU. For example, the relative mass growths of Saturn, Uranus and Neptune are somewhat 
ad-hoc, as well as their precise initial locations. Exploring the full range of possible giant 
planet histories and testing their effects on the inner solar system would be impossible. Instead 
we adopt the following strategy that we believe is more effective: we first investigate which 
ingredients of the scenario presented in the main paper have a significant role on the final results. 
We then neglect the aspects that do not matter and focus on those that do. 

More precisely, we will first address whether the presence of Uranus and Neptune has any 
role on the results. Then we address whether the assumed absence of type-I migration of Saturn 
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below 60 Earth masses has any role on the results as well as the actual chronology of its mass 
growth. 

Having concluded that the migration of Jupiter is the dominant aspect of the giant planet 
migration scenario that controls the sculpting of the asteroid belt, we will then explore in Sec- 
tion 5 how the results change with Jupiter's inward and outward migration rates and with its 
initial location. 

Do Uranus and Neptune play a role? The migration scenario presented in the main text (here- 
after the nominal simulation) included Uranus and Neptune growing from 5 M e , and only mi- 
grating outward when Jupiter and Saturn approached. The C-type material was partitioned 
between the giant planets and exterior to the orbit of Neptune. To understand whether Uranus 
and Neptune and their specific evolutions play a relevant role in the final results, we enact the 
same evolution as in the main paper for Jupiter and Saturn but discard Uranus and Neptune; the 
C-type material now resides in a disk exterior to the orbit of Saturn, from 6-10 AU. This range 
corresponds to the Uranus-Neptune belt and the trans-Neptunian disk in the nominal simulation. 

The S-type disk is scattered outward at a similar rate as in the nominal simulation, ~ 1 1%, 
and then ~ 0.7% is implanted into the asteroid belt during the outward migration of Jupiter and 
Saturn. The overall implantation efficiency is 0.87 x 10~ 3 , very similar to that found in the main 
text when Uranus and Neptune are included (0.7 x 1CT 3 ). From the C-type disk, ~ 0.5% of the 
planetesimals are scattered inward and are implanted in the asteroid belt. 

To compare, in the nominal simulation, the implantation efficiency was ~ 0.5% from the 
6-8 AU region (the Uranus -Neptune belt) and 0.02% from 8-13 AU. From this, we conclude 
that the presence of Uranus and Neptune just reduces the efficiencies slightly, particularly for 
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the most distant particles, which is reasonable as these particles have to undergo a more com- 
plex chain of planetary encounters before reaching the asteroid belt region. The final orbital 
distributions of the trapped asteroids are qualitatively similar in the nominal simulation and in 
the simulation without Uranus and Neptune (see Fig. S2). Based on these results from now on 
we exclude Uranus and Neptune from our tests, thus reducing the size of parameter space to 
investigate. 

Would type-I migration of Saturn play a role? Until recently, it was thought that small mass 
planets migrate inward (the so-called type-I migration 39 ) more rapidly than giant planets 40 . In 
this case, Saturn's core is expected to migrate towards Jupiter until it is captured in a mean mo- 
tion resonance with it; then Saturn's core follows the migration of Jupiter, staying in resonance 
while growing towards a Saturn-mass 54 . 

The scenario presented in the main paper incorporates modern results according to which 
rapid type-I inward migration does not occur in disks with a realistic cooling timescale until 
a mass of about 1/2 of a Saturn-mass is achieved 43 ' 45 ' 72 73 . However, radial migration is not 
necessarily obsolete. Planetary cores can migrate inward or outward at various rates depending 
on the properties of the disk (opacity, entropy gradient, viscosity; all possibly changing with 
time) and planetary masses. 

To test what effects the actual radial evolution of Saturn's core would have on the results, 
we adopt an extreme case: that type-I migration occurs as expected in isothermal disks. Thus, 
we test two scenarios inspired from hydrodynamical simulations 54 . 

- Saturn's core growing in the 1:2 resonance with Jupiter It has been shown that Saturn's 
core could be trapped in the 1:2 resonance with Jupiter, where it would accrete most of its 
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mass 54 . When Saturn is close to its final mass, the torque that it suffers from the disk becomes 
too strong relative to the resonant torque. At this time, the planet leaves the 1:2 resonance and 
very rapidly migrates into the 2:3 resonance. The capture in the 2:3 resonance reverses the 
migration direction of Jupiter and the two planets start to migrate outwards as in our nominal 
simulation. We implement this scenario as shown in Fig. S3. 

In this scenario, the inward migration of the planets scatters ~ 14% of all bodies beyond 
Saturn, and then ~ 0.8% of these are scattered back into the asteroid belt. The total efficiency is 
therefore 1 . 1 x 1CT 3 , very similar to scenario of Section 4 where Saturn's core had no migration. 

We do not need to test the capture of planetesimals from the outer disk in this scenario be- 
cause this disk is unaffected during the inward migration of the planets and, during the outward 
migration, this scenario and that of Section 4 are identical. 

- Saturn's core growing in the 2:3 resonance with Jupiter 

This scenario assumes that the core of Saturn was in the 2:3 resonance with Jupiter during the 
entire extent of its growth; it accreted gas until the planet became massive enough to force 
Jupiter to tack (Fig. S5). 

The simulation enacting this planetary evolution found that ~21 % of the bodies were 
scattered outward from the S-type disk during the inward migration, and subsequently ~1 .2% of 
those were scattered back into the asteroid belt region. Thus the entire efficiency was 2.5 x 1CT 3 , 
higher than in the previous cases, but of the same order of magnitude (all efficiencies for various 
particle sizes are reported below in Table S3). As before, we don't need to study the fate of the 
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outer disk because it would be identical to that discussed in Section 4. 

These results support that the actual radial evolution of Saturn's core does not significantly 
influence the final results. So, we assume, from now onwards, the simple scenario where Saturn 
is always in the 2:3 resonance with Jupiter, for simplicity. 

Chronology of Saturn's mass growth In the scenario of the main paper and in all tests above, 
the chronology of the mass growth was imposed on Saturn, and designed to create a "realistic" 
mass growth. However, given the vast uncertainties, this is essentially a free parameter. But, 
does the mass growth of Saturn change the results? 

To test this, we explore a scenario in which Saturn is in the 2:3 resonance with Jupiter 
and has its current mass from the beginning of the inward migration phase. We stress that this 
scenario has no physical meaning, because hydrodynamical simulations show that when Saturn 
is fully grown in the 2:3 resonance with Jupiter, the planets have to migrate outwards. However, 
this scenario is instrumental for our purpose, i.e. to test whether the actual mass history of 
Saturn influences the results or not. 

We also introduce a change in the migration history of Jupiter. Given that the inward 
migration might slow down before the outward migration starts, as anticipated in Sect. 2 we 
adopt a functional form for inward migration that is exponential rather than linear. For the 
inward migration we use r = 10 5 y, as for the outward migration. Jupiter switches to outward 
migration after 5 x 10 5 years from the beginning of the simulation. 

We find that the fraction of the S-disk objects ultimately trapped in the asteroid belt is 
0.8 x 1CT 3 , very similar to the trapping efficiencies found in all previous tests (~ 1 x 1CT 3 ). 
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The asteroid orbital distribution is also very similar to all other scenarios (Fig. S7). 

From all these tests we conclude that the basic feature in our model that determines the 
results discussed in the main paper is that Jupiter migrates inwards and then outwards, with the 
migration reversal when the planet is at 1.5 AU. The results of the asteroid belt do not depend 
significantly on how the cores of the other planets are assumed to grow or migrate. 

5 Exploring the parameter space further 

To have the reversal of Jupiter's migration it is essential that Saturn eventually enters in the 2:3 
resonance with Jupiter. We have discussed in Section 3 the conditions for this to happen, and 
we have argued that capture in resonance is required to explain the current structure of our giant 
planet system. 

There is no reason a priori that Saturn has to force Jupiter to reverse migration when 
Jupiter is at 1.5 AU. This could have happened at different distances. As explained in the main 
paper, the distance from the Sun at which Jupiter "tacked" needs to be determined by looking at 
constraints, such as the architecture of the terrestrial planet system and of the asteroid belt; the 
fact that the terrestrial planets are best reproduced when the disk of planetesimals from which 
they formed is truncated with an outer edge at 1 AU suggests that the tack occurred around 
1.5 AU. 

We now address the dependence of the results on two parameters that are expected to have 
a large influence: the initial location of Jupiter and its migration timescales. 

Given that the results are insensitive to the dynamics of the planetary cores and on the 
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chronology of Saturn's growth, for simplicity, we adopt the last model described above, where 
Saturn is always in the 2:3 resonance with Jupiter and has its current mass throughout the 
simulation. 

Initial location of Jupiter We present two additional simulations, one with Jupiter initially 
located at 2.5 AU and one at 4.5 AU. The inward and outward migration timescale r is still 
10 5 years. Only 100 km particles are considered. 

For the simulation with Jupiter at 2.5 AU, the S-type disk initially has an outer edge at 
2 AU and its total mass is 2.58 M e (in this case, the particle resolution was not decreased at 
all since they are all interior to 2.0 AU), with 14 embryos accounting for half of this mass. For 
the simulation with Jupiter at 4.5 AU, the disk initially extends out to 4.0 AU, with a total of 
4.73 M e , half of which is in 25 embryos. 

Table SI shows the capture probabilities for the S-type planetesimals. They increase 
monotonically with the initial location of Jupiter. This result was expected, in view of the ten- 
dency of planetesimals to be re-implanted near their original location at the end of the outward- 
then-inward scattering process. For the same reason, the initial location of Jupiter has a pro- 
found impact on the final semi major axis distribution of the trapped S-type planetesimals, as 
shown in Fig. S8. The distributions obtained with Jupiter starting at 2.5 or 3.5 AU give most 
of the objects inside of 2.8 AU (see also Fig. 3 of the main paper), in good agreement with 
observations (Fig. S9). The distribution with Jupiter starting at 4.5 AU is skewed toward large 
semi major axis. Thus, we tend to exclude that Jupiter formed so far out and we favor an initial 
location in the 2.5-3.5 interval, in good agreement with modern models of the location of the 
snow-line 70 ' 74 ' 75 . 
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Initial o, t,,^ 


± tot 


2.5 


2.81e-4 


3.5 


7.91e-4 


4.5 


1.37e-3 



Table SI: Probability that S-type planetesimals are captured in the asteroid belt, as a function 
of the initial location of Jupiter. The first column reports the initial semi major axis of the giant 
planet, in AU. The second column shows the final trapping probability P tot (to be compared 
with those reported in Table S3 and throughout the previous sections and main text). 

Again, we do not need to investigate the fate of planetesimals from the C-type disk in 
these simulations. The latter become active only during the outward planet migration, which is 
the same in these 3 simulations. 

Planet migration timescales Here we fix Jupiter's initial location at 3.5 AU, as in the refer- 
ence simulation of the main paper, and investigate the dependence of the results on the migration 
timescale r. Although we have done a broad exploration using different timescales for the in- 
ward and outward migration, for simplicity we present only the case where the two timescales 
are equal and focus on r = 10 5 , 2 x 10 5 and 5 x 10 5 y. In general, the results are more sensitive 
to the outward migration timescale. In fact, the speed at which Jupiter moves away from the 
asteroid belt is a crucial feature for the implantation process of both S-type planetesimals (pre- 
viously scattered outward during inward migration) and C-type planetesimals. Thus, presenting 
results for equal inward and outward migration timescales provides the significant insight into 
the range of outcomes. 

The probability that S-type planetesimals are implanted in the asteroid belt as a function of 
r is given in Table S2. As one sees, the capture probabilities for both planetesimal populations 
drop almost linearly with increasing migration timescale. This is expected because implantation 
of a planetesimal in the asteroid belt is more efficient if Jupiter moves away more rapidly after 
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having scattered the object inward. 



r 



P tot (S-type) P (C-type) 



le5 
2e5 
5e5 



7.91e-4 6.70e-3 
4.23e-4 4.37e-3 
1.98e-4 4.86e-3 



Table S2: Probability that S-type and C-type planetesimals are captured in the asteroid belt, as 
a function of the migration timescale r. The first column reports the r used for inward and 
outward migration. The second column shows the final trapping probability P tot for S-types 
(to be compared with P 2 reported in Table S3), and the third column reports the final trapping 
probability for the C-type disk between 6-9 AU (to be compared with the last three rows of 
Table S3, and 0.5% for the nominal simulation without Uranus and Neptune). 

From these results, the migration timescale remains quite unconstrained. Even a capture 
probability reduced by a factor of 5 would be sufficient, given the uncertainties on the amount of 
mass that the asteroid belt lost in its post-implantation evolution (see Section 7). Moreover, the 
original mass of the planetesimal disks is uncertain, particularly that of the outer disk. Thus, in 
absence of better constraints, we consider acceptable all the tested migration timescales. We re- 
mind the reader that hydrodynamical simulations show that r = 10 5 y is a reasonable timescale. 
But these timescales depend on disk's viscosity and scale height, so that hydrodynamical simu- 
lations cannot be used to predict the specific value of r for this problem. 

6 Trapping efficiencies as a function of planetesimal size 

All results illustrated thus far concern only planetesimals 100 km in size. Here we present 
the results for objects of 10 and 1000 km as tested in both cases where Saturn's core grows in 
resonance. Table S3 shows the probability for planetesimals originally in the S-type disk (inside 
the original location of Jupiter) to be trapped into the asteroid belt for the case with Saturn's 
core growing in the 2:3. Table S3 presents the same probabilities for planetesimals in the C- 
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type disk (between 6 and 9 AU). These tests can be compared directly to the results for those 
simulations done with 100 km particles. 



The final capture probabilities are weakly dependent on the particle sizes. The reason is 
that gas drag is not very important in the trapping process. This can be understood as follows: 
the planetesimals cannot be trapped permanently in the asteroid belt if Jupiter is still in the 
belt; so, they can be trapped only as Jupiter moves away, which occurs towards the end of the 
simulation. Since the gas decays exponentially with time during the outward migration of the 
giant planets there is not much gas left when the conditions are favorable for the trapping. 



Size 


Pi 


P 2 


Ptot 


P (C-type) 


10 km 


0.028 


0.0371 


1.04e-3 


0.0277 


100 km 


0.186 


0.0133 


2.47e-3 


6.7e-3 


1000 km 


0.191 


0.0050 


0.95e-3 


5.0e-3 



Table S3: Probability P tot to be captured in the asteroid belt for S-type and C-type planetesi- 
mals as a function of size. The first column reports the size of the particles. The fourth column 
Ptot is the product of the probability P 1 that an S-type object is scattered onto an orbit that is 
beyond that of Saturn and with a semi major axis within 9 AU and the probability P 2 that a 
particle in this S-type "scattered disk" is trapped into the asteroid belt during the planets' out- 
ward migration The fifth column is the probability for an object from the C-type disk between 
6-9 AU is implanted into the asteroid belt, which can be compared to P 2 . Percentages should 
not be affected by small number statistics: the probability Pi is computed on a population of 
4830 objects; then, the planetesimals trapped in the scattered disk are cloned 5-20 times, so that 
the probability P 2 is computed on a population of 2,000-7,000 objects. 



If gas drag is not the main mechanism, how are the planetesimals trapped in the asteroid 
belt? The main mechanism is that Jupiter ultimately scatters them inwards and leaves them be- 
hind when it migrates away. Most particles are observed to decouple dynamically from Jupiter 
when they are in the vicinity of its interior 2: 1 mean motion resonance. Interactions with mov- 
ing resonances can play a role in decoupling the planetesimals but these effects are, again, size 
independent. 
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Consistent with the explanation above, we have checked that the final orbital (a, e, i) dis- 
tributions of the trapped particles also depend weakly on size (Figure S10). 

7 Discussion on the Asteroid Belt 

The model presented in the main manuscript postulates asteroids being supplied to the main 
asteroid belt from two main parent regions, one from ~ 0.7-3.0 AU, and another from between 
or beyond the giant planets. Given that this runs counter to the classical idea that most or all 
of the asteroids were formed relatively close to where they currently reside, this demands some 
explanation of why our proposal is consistent with constraints. Thus, this is not meant to be a 
proof of two parent populations, rather a discussion of results which support this possibility. 

First and foremost, within the two populations represented in our model, the "S-types" and 
"C-types", we expect a compositional gradient to exist in each, as both populations originally 
extend radially over several AU. Thus these labels are used to simplify our model, not to suggest 
that the plethora of taxonomic types were all implanted from elsewhere. For example E- and S- 
types (respectively enstatite and ordinary chondrites) have much more in common than S- than 
C-type. The same is true for the C-, P- and D-types. We interpret E- and S-types as tracers of a 
compositional gradient in the volatile depleted disk inside of Jupiter's orbit, with E-types being 
originally closer to the Sun (maybe in the 1-2 AU zone), some of which find their way into the 
asteroid belt in our simulation. Note that isotopically, enstatite chondrites are the most similar 
to the Earth. Likewise, we interpret C-, P- and D-types as tracers of a compositional gradient 
in the giant planet region and beyond. Note that it has been proposed that P- and D-types are 
captured from the trans-Neptunian disk at the time of the late heavy bombardment 78 , i.e. from 
the outer edge of the primitive disk ~600 Myr after the S- and C-types are captured in this 
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scenario. These P- and D-type asteroids are the most similar to dormant comets, which makes 
sense if they are captured from the furthest reaches of the trans-neptunian disk. 

One could point to many differences between C-type and S-type asteroids as well as be- 
tween carbonaceous and ordinary chondrites (assuming that C-types and S-types are the parent 
bodies of carbonaceous and ordinary chondrites, respectively). From an astronomical point of 
view, C-types often show hydrated minerals in their spectrum, whereas S-types do not 79 . Some 
(Themis, and maybe Cybele) show water ice on their surface 80 ' 81 . When freshly cratered, some 
show cometary activity (e.g. the main belt comets) and when heated, some also show activity 
(e.g. Phaeton which is strictly a B-type, but similar to C-types). None of these features are 
common for S-types. Moreover, the bulk density of C-types (e.g. Mathilde) is about half that 
of S-types (e.g Gaspra, Ida, Eros) 82 , and the albedo of C-types is much lower than S-types. 

As for the meteoritic evidence, carbonaceous chondrites are much richer in water than the 
ordinary chondrites 83 . Also, hydration features are common among carbonaceous chondrites 
and not among the ordinary ones, which again shows a much larger presence of water in the 
parent bodies of the carbonaceous chondrites, at least in the past 84 . From a mineralogical point 
of view, in carbonaceous chondrites the matrix is much more abundant and CAIs are much 
more frequent, whereas type II chondrules are much more rare than in ordinary chondrites 85 . 
But perhaps the most significant difference between the two types of chondrites is their isotopic 
compositions. For both Chromium and Oxygen the carbonaceous chondrites have a continuum 
distribution in the isotopic diagram, whereas the ordinary chondrites are found in their own 
grouping 86 . Interestingly, as we do not find asteroids that are intermediate between C-type and 
S-type asteroids from a spectroscopic point of view, we also do not have meteorites that are 
intermediate between carbonaceous and ordinary chondrites from an elemental, mineralogical 
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and (particularly) isotopic point of view. 

These sharp differences can be explained if carbonaceous chondrites formed in a sig- 
nificantly different environment than the ordinary chondrites, so either in a different location 
of the proto-planetary disk or at a different time (assuming that the temperature and the chemi- 
cal/isotopic properties of the disk changes with both distance and time). In reality, though, there 
is no indication from any radioactive chronometer that the carbonaceous chondrites formed at a 
different time than the ordinary chondrites. In fact, new work that has measured the individual 
formation ages of chondrules in both carbonaceous and ordinary chondrites shows virtually the 
same age distributions 87 . So, if systematic differences in accretion times are ruled out, then 
it is necessary that the carbonaceous chondrites formed in a significantly different location of 
the disk, presumably further out from the Sun, where the temperature was lower (as evidenced 
by their richness in water). If the carbonaceous chondrites formed further out in the disk, the 
oxygen isotope variations would naturally be explained by the CO self-shielding process in the 
solar nebula 88 89 . 

The scientific picture of what comets really are is rapidly changing . The close flyby im- 
ages of comets (e.g. Comet Borrelly) show very little surface ice and small active regions 90 . 
The Stardust samples turned out to be quite similar to meteoritic samples 91 . Modeling work 
on the origin of the dust that produces the zodiacal light 92 predicts that at least 50% of the 
micro-meteorites collected on Earth are cometary; however, we see no clear separation of 
micro-meteorites into two categories, which could be traced to asteroidal and cometary dust 93 . 
Obviously, when we sample comets (and possibly also asteroids) through meteorites, micro- 
meteorites and aerogel traps, we sample only the rocky component, and not the ice. However, 
the continuum in physical properties between primitive asteroids and comets has recently been 
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well described 94 . 

Note however that we are not saying that C-type asteroids and comets are "identical". 
C-type asteroids in our models come from the 6-13 AU region. Current comets, from both 
the scattered disk and the Oort cloud, probably sample the outer extreme of the primitive disk, 
between 25-35 AU. There is still a significant distance between the source regions of C-type 
asteroids and comets, with many condensation fronts to cross (the CO front, for instance), so 
we still expect C-type asteroids to differ from comets in our proposed scenario. 

In principle, one could preserve the classical view of the asteroid belt within our scenario 
for the migration of the giant planets. This can be achieved with Jupiter (and Saturn) starting 
further out than in our nominal simulation. In fact, because bodies initially inside the orbit of 
Jupiter are re-captured in the asteroid belt with a statistical correlation between initial and final 
semi-major axes, we would find in the asteroid belt a compositional gradient if the disk inside 
of Jupiter (a) had such a gradient initially and (b) it extended originally out to ~3.5 AU. This 
would imply that Jupiter "formed" no closer than 4-4.5 AU. At the same time, the inner edge of 
the trans-Saturnian disk would be moved further out, implying few captured bodies from that 
part of the disk into the asteroid belt. This scenario is less appealing to explain the asteroid belt 
physical dichotomy, as it requires that all differences between the S-type and C-types arise over 
a relatively small semi major axis range, but it is in principle viable. 

8 What happened next? 

The solar system that emerges from the inward-then-outward migration of the giant planets is 
clearly not the solar system that is observed today. The giant planets are still in a compact, fully 
resonant configuration. The eccentricity distribution in the asteroid belt is skewed towards large 
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eccentricities. 

The evolution of the solar system after the disappearance of the gas disk has been investi- 
gated in detail over the last 5 years, and the scenario that we propose here is consistent with the 
understanding derived from those studies, reviewed below. 

Once the gas is removed, the giant planets still feel the gravitational interactions with the 
planetesimals that remain in a massive primitive disk outside of the orbit of Neptune. These 
interactions will eventually extract the giant planets from their mutual resonances. The planets 
then become unstable, and undergo a short phase during which they have mutual close en- 
counters. During this phase, Jupiter and Saturn acquire their current orbital eccentricities and 
inclinations. Uranus and Neptune are scattered outward on very eccentric orbits, which desta- 
bilize the distant planetesimal disk as a whole. The interaction with this dispersing disk damps 
the eccentricities of Uranus and Neptune by dynamical friction and eventually drive these plan- 
ets towards their current orbits. Recent iV-body simulations show that this kind of evolution is 
both possible and plausible 57 ' 70 . These works demonstrate that the giant planet instability is the 
dynamical link between the planets current configuration and the resonant orbital configuration 
that the giant planets must have had at the end of the gas-disk phase. 

When all these events happened is still a subject of debate. It is tempting to think that the 
instability of the giant planets occurred late, about 3.9 Gy ago or ~ 600 My after gas removal, 
because this would be a formidable explanation of the so-called Late Heavy Bombardment 
(LHB) of the terrestrial planets that is proposed to have occurred around that that time 95 . In- 
stead, if the giant planet instability occurred soon after the disappearance of the gas, some other 
explanation for the LHB needs to be found (see for instance the idea of a fifth planet stranded 
between Mars and the asteroid belt 96 ). Note that the very existence of the LHB is still debated 
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in the literature 97 . 

Regardless of whether the instability occurred early or late, its effects on the orbital distri- 
bution of the asteroid belt have already been investigated in great detail 69 98100 . To fix notation, 
we use the term "original distribution" to denote the orbital distribution that the asteroid belt 
had before the giant planet instability. Thus, if the model presented in this paper is correct, the 
"original distribution" would be similar to that shown in Fig. 3 of the main text. 

The works referenced above consistently demonstrate that Jupiter and Saturn had to evolve 
towards their current orbits very quickly, on a timescale shorter than 1 My, otherwise the orbital 
distribution of the asteroid belt would be very different from the current one, independent of 
the belt's "original distribution". The instability phase of the giant planets is compatible with 
these timescales for moving Jupiter and Saturn from their original 2:3 mean motion resonance 
to their current orbital configuration 69 . During the giant planets orbital re-configuration, the 
eccentricity distribution and, to a lesser extent, the inclination distribution in the asteroid belt 
are re-shuffled. The semi major axis distribution, however, is not directly affected. Thus, the 
final eccentricity distribution in our present model, illustrated for instance in Fig. 3 of the main 
paper, should not be compared to the current one, as it presumably evolved during the giant 
planet instability. In fact, discussions in the literature have argued that the asteroids' "original 
distribution" that best reproduces the current distribution is one skewed in eccentricity towards 
the stability boundary of the belt 101 . This is in agreement with what is shown in Fig. 3 of the 
main paper. Moreover, for the most rapid re-shuffling events the asteroids' "original distribu- 
tion" had to be roughly uniform in inclination in the 0-20° interval 69 . Again, Fig. 3 of the main 
paper fulfills this requirement. This is a non-trivial result, which has been difficult to obtain in 
planet formation simulations. For instance, a model 102 of the primordial sculpting of the aster- 
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oid belt resulted in an inclination distribution skewed towards large values, inconsistent with 
other findings 69 . 

The amount of material removed from the asteroid belt during the giant planet instability 
was estimated to be a factor of two 69 , when starting from a uniform eccentricity distribution. 
This factor could be 2-3 times higher if the eccentricity distribution is closer to Fig. 3 or to 
the distributions tested in some later planet migration scenarios 101 . Moreover, in the interim 
between the gas-disk disappearance and the onset of the giant planet instability (particularly 
if the latter occurred late), half of the population could be removed by dynamical diffusion, 
particularly given that most of the population is close to the stability boundary. This justifies 
the statement in the main paper that the asteroid population is "depleted by 50-90% during the 
LHB event". 

Finally, it has been shown that the asteroid belt population presumably declined by a 
factor of 3 since the time of the giant planet instability, by chaotic diffusion 103 . This process 
should have reduced the asteroid population without altering significantly its overall orbital 
distribution. 

In conclusion, we argue that the asteroid belt produced in the model presented in this paper 
is consistent with the current asteroid belt, given the uncertainties involved in its subsequent 
evolution. 
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Figure SI: The adopted gas profiles for three example distances of Jupiter: at 5.2 (solid), 1.5 
(dashed) and 3.5 AU (dotted). The amount of gas at any distance is then multiplied by an 
exponentially-decaying function of time, during the planets' outward migration. 
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Figure S2: Distribution of 100 km planetesimals at the end of giant planet migration for the 
scenario reported in Fig. 1 of the main text, but without Uranus and Neptune. The top panel 
shows the semi major axis distribution of the S-type bodies that will eventually end up in the 
asteroid belt, at three times: at the beginning of the simulation (dotted histogram), at the end of 
inward planet migration (dashed) and at the end of outward migration (solid). The second panel 
shows the final semi major axis distribution of the S-type (red histogram) and C-type (blue) 
planetesimals that are implanted in the asteroid belt. The bottom panels show orbital inclination 
and eccentricity vs. semi major axis. The dotted lines show the extent of the asteroid belt region 
for both inclination and eccentricity, and the dashed lines show the limits for perihelion less than 
1.0 (left line) and 1.5 (right line). Notice that the final distributions are similar to those shown 
in Fig. 3 of the main paper. 
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Figure S3: Top panel: the mass growth of Saturn in Earth masses, during inward migration. 
Bottom panel: evolution of the planets' semi major axes. Saturn is in the 1:2 mean motion 
resonance with Jupiter until it reaches its current mass. Then, it is extracted from this resonance 
and migrates in 10,000 years into the 2:3 resonance with Jupiter. This kind of evolution has 
been observed in hydrodynamical simulations 54 . 
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Figure S4: Distribution of 100 km planetesimals at the end of giant planet migration for the 
case of Saturn's core growing in the 1:2 resonance. The top panel shows the semi major axis 
distribution of the S-type bodies that will eventually end up in the asteroid belt, at three times: 
at the beginning of the simulation (dotted histogram), at the end of inward planet migration 
(dashed) and at the end of outward migration (solid). The second panel shows the final semi 
major axis distribution of the S-type (red histogram) and C-type (blue) planetesimals that are 
implanted in the asteroid belt. The bottom panels show orbital inclination and eccentricity vs. 
semi major axis. The dotted lines show the extent of the asteroid belt region for both inclination 
and eccentricity, and the dashed lines show the limits for perihelion less than 1.0 (left line) and 
1.5 (right line). Notice that the final distributions are similar to those shown in Fig. 3 of the 
main paper. 
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Figure S5: As in Fig. S3 but for the case where Saturn's core grows in the 2:3 resonance with 
Jupiter. 
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Figure S6: Distribution of 100 km planetesimals at the end of giant planet migration in the 
planetary evolution scenario of Fig. S5. The top panel shows the semi major axis distribution of 
the S-type bodies that will eventually end up in the asteroid belt, at three times: at the beginning 
of the simulation (dotted histogram), at the end of inward planet migration (dashed) and at the 
end of outward migration (solid). The second panel shows the final semi major axis distribution 
of the S-type (red histogram) and C-type ( blue) planetesimals that are implanted in the asteroid 
belt. The bottom panels show orbital inclination and eccentricity vs. semi major axis. The 
dotted lines show the extent of the asteroid belt region for both inclination and eccentricity, and 
the dashed lines show the limits for perihelion less than 1 .0 (left line) and 1 .5 (right line). Notice 
that the final distributions are similar to those shown in Fig. 3 of the main paper. 
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Figure S7: Distribution of 100 km planetesimals at the end of giant planet migration in the 
planetary evolution scenario with no mass growth for Saturn. The top panel shows the semi 
major axis distribution of the S-type bodies that will eventually end up in the asteroid belt, at 
three times: at the beginning of the simulation (dotted histogram), at the end of inward planet 
migration (dashed) and at the end of outward migration (solid). The second panel shows the final 
semi major axis distribution of the S-type (red histogram) and C-type ( blue) planetesimals that 
are implanted in the asteroid belt. The bottom panels show orbital inclination and eccentricity 
vs. semi major axis. The dotted lines show the extent of the asteroid belt region for both 
inclination and eccentricity, and the dashed lines show the limits for perihelion less than 1 .0 
(left line) and 1.5 (right line). Notice that the final distributions are similar to those shown in 
Fig. 3 of the main paper. 
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Figure S8: Final distributions of S-type planetesimals at the end of the inward-then-outward 
migration process. Different colors denote the initial location of Jupiter in the corresponding 
simulation: 2.5 AU (red), 3.5 AU (black) and 4.5 AU (blue). Top panel: the distribution on 
the eccentricity vs. semi major axis plane. Filled circles correspond to particles in the asteroid 
belt region at the end of the simulation (with inclination below 25°). The only particles plotted 
are those which are scattered back inwards from the 3-9 AU region (hence the planetesimals 
and embryos between 0.5-1.0 are not shown). Bottom panel: normalized histograms for the 
particles which are placed in the asteroid belt region. 
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Figure S9: Number of observed asteroids larger than 100 km in diameter that are C-type (blue 
histogram) and S-type (red histogram), binned as a function of semi major axis 76,77 
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Figure S10: Orbital distribution of particles for the scenario of Saturn growing in the 2:3 res- 
onance with Jupiter for 10 km and 1000 km bodies. The top two panels show the eccentricity 
and inclination as a function of semi major axis for 1000 km bodies. The bottom two panels 
show the same properties for 10 km bodies. In all panels, bodies within the asteroid belt region 
are filled symbols, and those outside are open symbols. 
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